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The search for Majorana fermions in p-wave paired fermionic systems has recently moved to the 
forefront of condensed-matter research. Here we propose an alternative route and show theoreti- 
cally that Majorana-like modes can be realized and probed in a driven-dissipative system of strongly 
correlated photons consisting of a chain of tunnel-coupled cavities, where p-wave pairing effectively 
arises from the interplay between strong on-site interactions and two-photon parametric driving. 
The nonlocal nature of these exotic modes could be demonstrated through cross-correlation mea- 
surements carried out at the ends of the chain — revealing a strong photon bunching signature — and 
their non-Abelian properties could be simulated through tunnel-braid operations. 

PACS numbers: 42.50.Pq, 03.75.Lm, 73.21.-b 



In recent years, strongly correlated photons have 
proved to be a remarkably rich platform for investigating 
phenomena traditionally regarded as pertaining to con- 
densed matter physics. Tremendous theoretical and ex- 
perimental efforts have made it possible to achieve strong 
optical nonlinearities at the single-photon level [1, 2] and 
to demonstrate photon blockade effects [3-8]. Meanwhile, 
the pursuit of Majorana fermions has become a new focus 
of condensed matter research [9, 10], and p-wave paired 
superfluids and superconductors have been promoted as 
paradigmatic systems for investigating the physics of Ma- 
jorana modes [11]. Even though strongly interacting pho- 
tons have been predicted to exhibit a typical fermionic 
behavior [12, 13], optical systems have never been con- 
sidered as candidates for realizing such exotic physics. 

In this Letter, we show that Majorana-like modes 
(MLMs) can be obtained in a one-dimensional (ID) 
strongly correlated system of impenetrable (or "fermion- 
ized") photons. More specifically, we consider a chain 
of coupled cavities with strong on-site nonlinearities and 
introduce a drive mechanism based on parametric am- 
plification which, in stark contrast to previous works, 
gives rise to an effective p-wave pairing between (fermion- 
ized) photons. We map our system to the ID chain 
originally proposed by Kitaev as a toy model for Ma- 
jorana fermions [14], and show the existence of zero- 
energy modes with properties similar to those of Majo- 
rana modes in solid-state systems. Owing to the intrin- 
sic dissipative nature of the system, these "Majorana- 
like" modes do not benefit from parity (or "topological" ) 
protection against decoherence [15, 16], and thus can- 
not serve as topological quantum memories [17]. How- 
ever, they do behave as genuine Majorana modes on time 
scales shorter than the lifetime of a photon in the system, 
allowing Majorana physics to be probed. 

To demonstrate the fact that MLMs can be detected 
via simple optical schemes, we propose a realistic exper- 
iment that takes full advantage of the optical nature of 
the system and allows for the direct observation of MLMs 
through second-order photon cross-correlation measure- 
ments. Although our proposal is strictly limited to ID — 



since impenetrable photons do not behave as fermions in 
higher dimensions — we show that MLMs can effectively 
be exchanged using "tunnel-braid" operations [18], en- 
abling us to simulate their non-Abelian properties. 

The model. — The backbone of our system consists of a 
ID chain of N optical cavities coupled through nearest- 
neighbor photon tunneling [Fig. 1] . Each cavity exhibits 
a large optical nonlinearity (i.e., is strongly coupled to 
an artificial atom) and supports a single mode that can 
be described as a Wannier function localized on site i 
around the cavity center. Photon tunneling occurs as 
a result of the non-vanishing spatial overlap between 
nearest-neighboring Wannier modes [19], and the system 
Hamiltonian takes the generic form of a generalized Bose- 
Hubbard model: 

N N N-1 

= 5Z bjb, + I ^ blblbA -JJ2 (^l^'+i + ^■^•)' 

2 — 1 i—1 2—1 

(1) 

where bi (b^) are annihilation (creation) operators asso- 
ciated with the i"^ cavity of the chain with resonance 
frequency is the strength of the on-site photon- 

photon repulsion (Kerr energy) due to the large optical 
nonlinearities, and J denotes the tunneling amplitude be- 




FIG. 1. (Color online). Driven-dissipative chain of coupled 
cavities. Each cavity exhibits a large optical nonlinearity and 
sustains a single Wannier mode which weakly overlaps with 
its nearest neighbors, thus allowing for photon hopping be- 
tween sites. Parametric pumps (depicted by arrows) couple 
to the weak inter-cavity field and inject photon pairs which, 
owing to strong photon-photon repulsion, split up into differ- 
ent cavities, effectively giving rise to p-wave pairing. 



2 



twccn nearest-neighboring sites. In this work, we will fo- 
cus exclusively on the strong-interaction regime, in which 
the energy cost U of adding an extra photon to an oc- 
cupied cavity is by far the largest of all relevant energy 
scales [20]. In this so-called "hard-core" limit, the oc- 
cupation of each site is effectively restricted to or 1, 
and the photons exhibit a characteristic fermionic be- 
havior [12, 13]. 

To achieve p-wave pairing, we introduce parametric 
pumps (or amplifiers) which, in stark contrast to usual 
coherent drives, are tailored to inject pairs of photons 
into the system through nonlinear optical processes (see 
Supplemental Material [21]). Assuming that these pumps 
drive the system locally through the inter-cavity field — 
which consists of a superposition of two neighboring 
Wannier modes — photons from a single pair can either 
be emitted into the same cavity (or Wannier mode), or 
settle into different, nearest-neighboring cavities. In the 
strong-interaction regime, the second process is strongly 
favored, and the drive Hamiltonian effectively reads 

N-l 

/fdrive = -I A] 5^ {e'^^^-^+^^bA+i + h.c), (2) 

i=l 

where A = )A)e"^ defines the amplitude and phase of the 
parametric pumps, and ojp their frequency. We note that 
the amplitude |A] of the parametric drive is determined 
by the overlap of the Wannier modes in a similar way as 
the tunneling amplitude J defined above. We thus expect 
to be able to reach a regime in which the two quantities 
are of the same order. Physically, the above Hamiltonian 
describes the coherent exchange ofp-wave paired photons 
between the system and the classical pump field(s) [22]. 
It provides the optical counterpart of p-wave supercon- 
ductivity that is crucial to access Majorana physics, and 
compensates for losses by continuously replenishing the 
system with photons. The time evolution of the system 
including the drive and the photon losses is governed by 
the Lindblad master equation 

N 

dtp = -i [Ho + J?drive,p] + (bipbl - l{blbi,p}') , 

i=l 

(3) 

where p is the density matrix of the system and T the 
decay rate associated with the individual cavities. 

Mapping to a ID Kitaev chain. — In the strong- 
interaction regime (U S> J, |A|), the Hilbcrt space of the 
system effectively reduces to that of hard-core photons 
bi = VbiP, b\ = Vb\V, where V projects onto the sub- 
space of single occupancy. Hard-core photons can be seen 
as spin-1/2 particles, with Pauli-type matrices a~ = 2bi, 
u'l = 2b\ [af = af±iaf), and their fermionic nature can 
be unveiled by mapping the spin-1/2 particles to spin- 
less fermions a^, aj using a Jordan- Wigner transforma- 
tion [23] of the form = 5 11^=1 ("'^l)^^- Defining 



p = ujp — iOc and moving to a rotating frame defined by 
Hi = WpX^j^i^i, the Hamiltonian H = Hq + iJdrive of 
the full system becomes, in the fermionic picture, 

iV-l N-l 

H = -JY^ {alai+i + h.c.) + \A\ ^ (e^'^aiOi+i -|- h.c.) 

i=l i=l 

N 

- p'^a\ai, (4) 

i=l 

which corresponds to the ID p-wave superconductor of 
spinless fermions originally introduced by Kitaev [14]. 
Here, the cavity frequency lOc plays the role of a Fermi 
level, and the detuning p = u)p — tOc between the pump 
and cavity frequencies that of a chemical potential. As- 
suming that I A] 7^ 0, two topologically distinct (gapped) 
phases can be identified [14]: a trivial phase correspond- 
ing to \p\ > \2J\, and a nontrivial phase corresponding 
to \p\ < \2J\, in which the system supports Majorana 
modes that are exponentially localized at both ends of 
the chain. The topological phenomena associated with 
the ID Kitaev chain can be most easily understood for 
J = A > and p = (i.e. uip = Wg). In this illustrative 
case, the Hamiltonian reduces to 

N-\ N-l 

i=l i-1 

with Majorana operators defined as 

C2i-i = ai + a] = n}=i(-<^j)<^?, (.g^ 
C2i =-i{ai-al) = -nj=l(-<^|)f^]'> 

and can readily be diagonalized as H = 2J^^"^^(ata^ — 
1/2) with Bogoliubov-Valatin quasiparticle operators 
fli = (c2i + ic2i-i-i)/2. The associated spectrum is sym- 
metric about the "Fermi level" ojc and features a gap 
2J. Most importantly, it exhibits two Majorana zero- 
energy modes corresponding to the Majorana operators 
ci and C2N localized at the ends of the chain but ab- 
sent from the Hamiltonian. These modes define a two- 
dimensional, nonlocal degenerate (zero-energy) subspace 
which we identify as a "Majorana qubit", with associated 
operator 

N 

(t'm = iciC2JV = l[{-(T])a^a^^. (7) 

The string-like operator P = njLi(~cr|) that appears in 
the above expression corresponds to the parity operator 
associated with the total number of (fermionized) pho- 
tons. It commutes with the Hamiltonian of Eq. (4), but 
anticommutes with the collapse operators entering 
the Liouvillian in Eq. (3). Physically, this means that 
single-photon losses result in the breakdown of parity 
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conservation, as expected, such that the Majorana qubit 
is not parity-protected — or "topologically protected" — 
from decoherence [24]. Although this disqualifies Ma- 
jorana modes of light for practical applications such as 
topological quantum memories, this crucially does not 
hinder the observation of their exotic physics, since the 
key features of A^lajorana physics — such as the existence 
of localized Majorana modes and the possibility of sim- 
ulating non-Abelian braiding operations in the associ- 
ated zero-energy subspace — do not require perfect parity 
protection to be realized. In fact, as we demonstrate be- 
low, Majorana physics does remain accessible within time 
scales much shorter than the lifetime ^ l/F of a photon 
in the system. However, the Majorana modes appearing 
in our optical framework only behave as genuine, local 
Majorana modes in the limit F — > 0. Only then does 
the string-like operator P that they carry — a remnant of 
the Jordan- Wigner mapping — reduce to a simple phase 
P = ±1. In this respect, we will refer to Majorana modes 
of light as "Majorana- like" modes (MLMs). 

Our proposal for MLMs, embodied in Eqs. (1), (2), 
and (3) , could be realized in any cavity quantum electro- 
dynamics (cavity QED) system lying deep in the strong- 
coupling regime. In the Supplemental Material [21], we 
outline a potential implementation in circuit QED which 
we deem as closest to experimental realization, and give 
ballpark figures for the relevant parameters and energy 
scales. 

Optical detection scheme. — Multiple schemes have 
been proposed for detecting Majorana modes in solid- 
state systems (sec, e.g., [10] for a review). Although we 
believe that most of them can be transposed to our op- 
tical setting, we will focus on the detection of Majorana- 
mediated (photonic) Cooper pair splitting, following the 
proposals of Refs. [25, 26]. We start with a Kitaev chain 
in a topologically nontrivial regime defined by J, A > 0, 
and < /i < 2J, without loss of generality. In such a 
parameter range, exponentially localized MLMs are ex- 
pected on both sides of the chain, with a length scale 
that increases with /i and diverges a.s fi ^ 2J [14]. For 
finite iJL <2J and small enough system sizes, these modes 
weakly couple and the levels of the Majorana qubit that 
they form split in energy by an amount 5m > [27], as 
captured by the Hamiltonian (5m (the explicit form 
of f7|^ reduces to cr^tr^ when restricted to the end cav- 
ities, and is given by Eq. (7) for /i = 0); as a key in- 
gredient, we assume that 6m ^ Eg, where Eg denotes 
the gap of the system [28]. Next we introduce two addi- 
tional nonlinear cavities — one on each side of the Kitaev 
chain — which we refer to as the left (L) and right (R) 
probe cavities, respectively. We assume that the latter 
have resonance frequencies WL,i?. = s^nd that they 
couple to the end cavities of the Kitaev chain through 
weak tunneling only. In the rotating frame introduced 
above (in deriving Eq. (4)), the Hamiltonian describ- 
ing the interaction with the probes then takes the form 



^fprobo = -JLi(7l(^f + (^lf7f)-JRi(T%(T'^ + alfal), where 
< Jl,r *C <5m denotes the weak amplitude for tunneling 
into the left and right probe cavities, respectively. Since 
all energy scales associated with -ffprobo are, by assiimp- 
tion, much smaller than Eg, the probe cavities only probe 
the low-energy physics associated with the MLMs of the 
chain. Owing to this energy selectivity, the terms of the 
form cfaj appearing in iJprobe — which mediate coupling 
to higher excited states — can safely be neglected, while 
the operators af and aff — which anticommute with (t|| 
and thus effectively describe a spin fiip of the Majorana 
qubit — can be replaced by a^. This results in the follow- 
ing low-energy effective Hamiltonian for the full system: 



(8) 



Physically, the above expression tells us that the non- 
local Majorana qubit formed by the "localized" (in the 
limit F — >■ 0) MLMs of the chain mediates a nonlocal 
coherent exchange of photons between the probe cavi- 
ties. Clearly, the bottleneck of such an exchange is given 
by the time scale tM ~ over which the Majorana 

qubit evolves. We thus only expect to see nonlocal cor- 
relations between the probes if Im is the shortest time 
scale in Eq. (8), i.e., if Sm '> Jl,r, as has been assumed. 
To detect these correlations, one can take advantage of 
the intrinsic dissipative nature of the system. Assum- 
ing that the decay rate of the probe cavities satisfies 
Tl./j ^ Jl.r ^ Sm, so that spontaneous emission oc- 
curs on a time scale much longer than the time scale tM 
over which correlations are generated, we expect a direct 
signature of MLMs to appear in the second-order pho- 
ton cross-correlations between the light emitted from the 
two probe cavities. In order to illustrate this, we consider 
the simple case Jl = Jr = \f2J. F^ = Yr = SF, in the 
limit where the decay rate F associated with the cavities 
of the Kitaev chain vanishes. Following the method of 
Ref. [29] (see Supplemental Material [21]), we then ob- 
tain steady-state photon cross-correlations between the 
probe cavities that read 



(2) ^ (^L^^i^^i.) 



1 + 



\" L iV'Rl 



(1 + (c7£))(l + (a^)) 



1 + 



(J2+F2)2 



(9) 



Remembering that f ^ ^ ^ ^2^^ \\\w& find g^^^ :g> 2; 
in other words, the light emitted by the spatially sepa- 
rated probe cavities is strongly bunched [30]. To exam- 
ine the effect of weak dissipation from the chain, we have 
carried out numerical simulations of the jull Kitaev chain 
coupled to the probe cavities. Our studies confirm that 
a striking nonlocal photon-bunching signature of MLMs 
remains visible for small decay rates F ~ Vj^^r <C 5m of 
the chain cavities and small enough system sizes, i.e., as 
long as the effective width of the Majorana levels is much 
smaller than their energy splitting [Fig. 2]. 
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FIG. 2. (Color online). Numerical results showing the second- 
order cross-correlation function g)^^ as a function of the de- 
tuning /i for A'^ = 12 (including probe cavities), A/ J — 1, 
Jl,r/J = 0.02, and for different values of F/J = Tl,r/J- 
MLMs (with energy splitting 5m determined by fi) are ex- 
pected in the region < jj,/ J < 2 (delimited by a dashed 
vertical line) which corresponds to a topologically nontrivial 
phase (in the limit A*' — > oo). As expected, a strong bunch- 
ing signature is observed for large /i inside the topologically 
nontrivial region, while g)^^ ~ 1 beyond the critical point 
fj,/ J = 2, clearly signaling the absence of MLMs. All results 
were obtained using a Monte Carlo wave-function approach 
(see, e.g., [31]) with 400 trajectories. 



We remark that the above results closely parallel those 
obtained in Refs. [25, 26] in the solid-state setting. Here, 
the probe cavities play a similar role as the metallic leads 
used in typical solid-state proposals: they provide a nar- 
row band of states of effective width ^ l/r^^fl close to 
the "Fermi level" ujc into which fermionized photons can 
be emitted from the system. Due to the large splitting 6m 
between MLMs, nonlocal — or "crossed" — Andreev reflec- 
tion is favored over its local analog, and Cooper pairs of 
photons are split into separate leads (or probe cavities). 
We remark that one could in principle remove the probe 
cavities and observe directly the light emitted by the end 
cavities of the Kitaev chain. In that case, however, spec- 
tral filtering would be required in order to isolate the 
bunching signature associated with the low-energy MLMs 
from contributions of the bulk. 

Versatility of the optical proposal beyond detection. — 
In contrast to previous proposals (such as Refs. [32-34]), 
our optical system provides a very versatile platform for 
simulating Majorana physics. In addition to being a con- 
ceptually simple realization of Kitaev's original lattice 
toy model — fermionized photons are intrinsically spin- 
less, and pairing occurs between nearest-neighboring cav- 
ities only — it allows for single-site addressability and for 
local control of all parameters entering Kitaev's model: 
the "chemical potential" /i can be tuned easily and locally 
by modifying the frequency of the individual pumps (or, 
alternatively, the resonance frequency of the individual 



cavities) , and the amplitude and phase of the "supercon- 
ducting order parameter" A can similarly be adjusted 
by controlling the amplitude and phase of the paramet- 
ric amplifiers. The tunneling amplitude J, on the other 
hand, can be tuned by introducing intermediate control 
devices between cavities (see, e.g., Ref. [35] and refs. 
therein). Such level of control is key to overcoming the 
crucial challenges currently facing most solid-state pro- 
posals, such as the tuning in and out of the topological 
phase and the suppression of disorder effects [34]. 

Despite its conceptual simplicity, versatility, and phys- 
ical realizability, our optical proposal departs from be- 
ing ideal in two respects: (i) it lacks topological pro- 
tection (or parity conservation), and (ii) it cannot be 
scaled up to networks of ID wires [36, 37]. The first 
imperfection arises as a direct consequence of photon 
losses — unavoidable in photonic systems — and puts strin- 
gent constraints on the time scale over which Majorana 
physics can be observed; namely, MLMs must be manip- 
ulated and detected on a time scale much shorter than 
the lifetime of a photon in the system. We argue in the 
Supplemental Material [21], however, that state-of-the- 
art technologies in circuit QED could allow for the ex- 
perimental realization of our proposal with a sufficient 
control over dissipation to meet these requirements. The 
second imperfection stems from the intrinsic nonlocal na- 
ture of the Jordan- Wigner mapping invoked in deriving 
Eq. (4) , rufing out the possibility to observe non- Abefian 
exchange statistics in connected wire geometries. In ID, 
the Jordan- Wigner string carried by the end MLMs of 
a chain essentially corresponds to the parity operator of 
the latter (see, e.g., Eq. (7)), such that end MLMs do 
behave as genuine, local Majorana modes on time scales 
over which parity is effectively constant. In higher dimen- 
sions, however, the situation changes drastically: when 
multiple ID chains are contacted (not through their ends, 
so that the systems effectively is higher-dimensional), 
the parity of the individual chains becomes a dynami- 
cal quantity and the nonlocal nature of the MLMs comes 
into play — with dramatic consequences such as the ab- 
sence of non-Abelian exchange statistics. To avoid such 
complications, we simply strictly restrict ourselves to ID 
systems. In that case, the exchange of Majorana modes 
is impossible in real space, but can nevertheless be sim- 
ulated using so-called "tunnel-braid" operations [18] . As 
shown in the Supplemental Material [21], these oper- 
ations only preserve the degeneracy of the MLMs — as 
real-space braiding — provided that the relative phase be- 
tween the latter is properly tuned, and therefore are not 
strictly speaking topologically protected. This, however, 
does not constitute an additional problem in our opti- 
cal setting where parity is anyway not conserved. In the 
framework of our proposal, tunnel-braid operations cru- 
cially allow us to obviate the need for real-space braiding, 
hence providing us with a full-fiedged ID optical platform 
for Majorana physics. 
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Conclusions. — Wc anticipate that our proposal for re- 
alizing and detecting photonic p-wave pairing will allow 
for an exciting alternative avenue for the investigation of 
Majorana physics. An interesting possibility would be 
the investigation of Majorana modes in a continuum ID 
model of strongly interacting optical photons [12]. 
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P-WAVE PAIRING OF STRONGLY CORRELATED PHOTONS 



In this section, we demonstrate explicitly that p-wave pairing generally emerges as a 
result of the interplay between two-photon parametric pumping and strong on-site photon- 
photon repulsion, and derive the effective drive Hamiltonian of Eq. (2) in the main text. 
Although multiple schemes can be envisioned to realize the key nonlinear pumping process 
that takes photons from a pump field and injects coherent photon pairs into the system, we 
will illustrate the general mechanism that leads to p-wave pairing in a model based on a 
theoretical description of coupled cavity arrays that was developed in Ref . [1] . 

We start by assuming that the parametric pumps act locally in between the nearest- 
neighboring cavities corresponding to the sites i and i + 1 according to the standard Hamil- 
tonian [2] 

i^drive. = / d'rx^'\r)Elf{r,t)E^^\r)E^^\r) + h.c, (1) 

where E^'^\r, t) is the positive-frequency part of the pump (index "p") optical field, £'g~^(r) 
the negative-frequency part of the generated signal (index "s" ) optical fields, and x^^^ (r) the 
effective second-order optical nonlincarity of the system (taking into account the polarization 
vectors of the three interacting fields). We assume that the pump field can be treated as 
a classical monochromatic field of frequency 2ujp^i and complex amplitude E^^-, and that 
there is no significant depletion of the latter such that E^ - can be considered as constant. 
Denoting the spatial mode function of the pump as (pp^i{r), Ep'^\r, t) then takes the form 

Ei-!,\r,t) = El^p,ir)e-''-^K (2) 

The generated signal optical fields, on the other hand, can generally be expanded in terms 
of the Wannicr modes 0i(r) of the coupled cavity array. Assuming that these modes decay 
sufficiently fast so that only nearest-neighboring modes have a non-vanishing overlap, we 
can write 

(r) = -i (0*(r)6l + 0*+i(r)6l+i) , (3) 
where ouc is the resonance frequency of the cavity (assumed to be the same for all cavities). 
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Plugging Eqs. (2) and (3) into Eq. (1), we thus obtain 

where we have invoked the presumably strong on-site repulsion between photons to neglect 
all terms of the form The strength of the degenerate parametric pumping therefore 
clearly appears to depend on the overlap of the interacting modes, as expected. Assuming, 
for simplicity, that the pump field and the effective second-order optical nonlinearity are 
essentially constant over the volume between the cavities, and that the pump field vanishes 
outside this volume, we finally obtain 

H^riye,i ~ (^^ dV 0*(r)0*^i (r)^ blbl.e-''^^' + h.c. 

^ -a; blbl.e-'-'^^' + h.c, (5) 

where the complex coupling constant /tj essentially depends (linearly) on x*^^) and on the 
amplitude of the pump field E^ -. This last result shows that the strength of the two- photon 
parametric drive (i.e., the p-wave pairing strength) directly depends on the overlap of the 
Wannier modes corresponding to nearest-neighboring cavities, as argued in the main text. 



STEADY-STATE SECOND-ORDER PHOTON CROSS-CORRELATIONS 

In what follows, we give a brief derivation of Eq. (11) from the main text. In accordance 
with the discussion presented there, we start with the Lindblad master equation 

dtp = -i [H,s, p]+Yl [^^P - P}) ' (6) 

i=L,R 

where the Hamiltonian and the Lindblad operators are respectively given by 

-f^eflf = Smctm - V2J {ala^ + ctmCt"^) , (7) 

Li = \ll<y-. (8) 

Regarding the system as a chain of 3 interacting spins and relabeling the sites as L — > 1, 
M — > 2, — > 3, one can perform a Jordan- Wigner transformation and describe the problem 
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in a basis of Majorana operators defined as 

C2i-i = (n 



j=i (^jK, 

=(n;;W)'^r- 



(9) 



In tfiis fermionic description, tfie Hamiltonian and tfie Lindblad operators take tfie form 

i^eff = -i5MC3C4 - i\/2J(c2C3 + C4C5), (10) 

L^^^^{c^-ic^), (11) 

^3 = -\/f(c5-iC6), (12) 

and can conveniently be expressed as H^q — ^c^Hc and Lj = Ijc, where — (ci, C2, . . . , Cg). 
Note that we have dropped the string-hke operator of the form 11^=1 ^^^^ would normally 
appear in L3, since the latter does not affect Eq. (6) (see, e.g., Ref. [3]). The fact that 
Eq. (6) is quadratic allows for a great simplification of the problem. Assuming, without 
loss of generality, that the initial state of the system has a Gaussian form, one can focus 
exclusively on the time evolution of the correlation matrix Cij — (i/2) tr([cj, Cj]p) which 
contains all information about the system. As shown in Ref. [4], the master equation of 
Eq. (6) then reduces to a matrix equation of the form 

dtC = X^^C + CX-Y, (13) 

where X = -2iH + 2 Re M and y = 4 Im M, with M = 1? ® ij. Here we find 

/fO oo\ 
f 2V2J 

-2V2J -5m 
5m 

\ 

and the steady-state correlation matrix satisfies 

X^C + CX^Y. (15) 

In order to solve the above matrix equation, we follow the method of Ref. [3] and make the 
ansatz 



X 



2V2J 

-2^27 f 

or/ 



Y 



/ 





2f 


















-2f 
















































































2f 




\ 














-2f 





/ 



(14) 



X 



T 



(16) 
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where and Wj are left and right eigenvectors of X, respectively. Introducing this ansatz 
into Eq. (15) and assuming that Vj, Wj are normalized according to v]vfe = Sj^k, we find 
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(17) 



Xj + Afc 

where Xj are the eigenvalues corresponding to the right (or left) eigenvectors of X. Therefore, 
solving for the steady-state correlation matrix basically amounts to solving two eigenvalue 
problems: = XjWj and VjX — XjVj. Using the above exphcit form of X and Y, 

straightforward algebra leads to 

/ n l^-r^+t^+l) %/27t(7^+t^) 727^t 7t' ^ 





7"(7"+t"+l) 

_V27t (72+42) 
" 72 + (72+t2)2 
V272f 
72 + (72+t2)2 

72 + (72+t2)2 





7^(7^+«^+l) 


7274(72 +t2) 


V2'yH 


7t^ 


72+(72+t2)^ 


72 + (72+t2)^ 


72+(72+t2)2 


72 + (72+t2)2 


















































7*2 




v^7i (72+42) 


7^(7^+4^+1) 


72 + (72+«2)2 


72+(72+*2)2 


72 + (72+t2)2 


72 + (72+t2)2 



7^2 

72 + (72+t2)2 

72724 

72 + (72+t2)2 
x/274(72 +42) 

72 + (72+t2)2 

7^(7^+*^+!) 

72+(72+t2)2 





(18) 

where we have defined the dimensionless quantities t — J /5m, J — i'/Su- The steady-state 
second-order photon cross-correlations between the left and right probes then read 



(2) 1 



= 1 + 
= 1 - 

= 1 + 



(l + (ai))(l + (a|j)) 

{CiC5){c2Cq) - {ciCq){c2C5) 
(1 - i(ciC2))(l - i(c5C6)) 
Cl,5C'2,6 — C'l,6C'2,5 

(l-Ci,2)(l-a,6) 

7^ 



(19) 



(^2+^2)2' 

where we have used Wick's theorem in the second equality. Similarly, the steady-state 
occupation of the probe cavities is given by 



1 + {^Ir) 



7' + {t^ + 7')' 



(20) 



(2) 

and scales as l/gL^ when J ~ F. 
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IMPLEMENTATION IN CIRCUIT QED 



Our proposal for accessing Majorana physics in the optical context requires two crucial 
ingredients: strong coupling and (comparatively) weak photon losses, both of which are cur- 
rently available in superconducting-circuit-based cavity quantum electrodynamics (circuit 
QED) (see Ref. [5], e.g., for a recent review). Superconducting qubits thus stand apart as 
an ideal platform for realizing our proposal — even more so because of the great flexibility in 
fabrication and control provided by state-of-the-art technologies. We outhne below a poten- 
tial implementation of our proposal which we deem as closest to experimental realization, 
and give ballpark figures for the relevant parameters and energy scales. 

The system that we typically envision consists of a chain of capacitively-coupled identical 
microwave resonators (playing the role of cavities) , each incorporating two superconducting 
transmon qubits (playing the role of artificial atoms) . The first of these qubits is positioned at 
an antinode of the intracavity field, with a cavity-qubit coupling strength g ~ oj^sf^i where 
lO}) is the bare-cavity resonance frequency and a the fine-structure constant; the strength of 
the effective on-site interaction — or, equivalently, of the anharmonicity — associated with the 
lower polariton mode (with frequency o'c = cjb — g) is then given by U = {2- ^p^)g. The 
second qubit is located in between pairs of neighboring cavities, and is driven parametrically 
using strong microwave fields at frequency luj^ so as to ensure the generation or annihilation 
of photon pairs at frequency oub — g. In the limit where the amplitude of the parametric 
drive (or "pump") is weak as compared to g, the only energetically allowed process is the 
generation (annihilation) of photon pairs in neighboring cavities, as discussed in the main 
text and in the first section above. This results in an effective p-wave pairing interaction 
with a phase determined by that of the microwave drive fields. 

We emphasize that g ~ O-luj, ~ lO^F, where F is the cavity decay rate (see main text), is 
routinely obtained in state-of-the-art circuit-QED devices. This implies that the tunneling 
amplitude J between neighboring cavities and the amplitude |A| of the parametric pumps 
can in principle be tailored to satisfy the relation U ~5'^J, |A|^r which lies at the core 
of our proposal. The first inequality {U ^ g ':$> J,\A\) ensures that the system lies deep in 
the strong-coupling regime where photons effectively behave as fermions (or "fermionized" 
photons) , and constitutes a sufficient condition for our model to map to that of Kitaev (see 
main text). If one additionally adjust the detuning /i — oup — ouc between the cavities and 
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the pumps so that < 2 J, the system is found in a topologically nontrivial regime and 
exhibits Majorana-hke modes with an energy sphtting 5m that depends on J, |A|, and |//|, 
and vanishes exponentially with the system size N (see the explicit form below) . The second 
inequality (J, |A| ^ F), in contrast, determines the broadening of the Majorana levels, i.e., 
how well they can be resolved for a fixed energy splitting Sm- For Majorana-like modes to 
be detectable using the scheme presented in the main text, this broadening must crucially 
be (i) much smaller than the energy splitting 5m between the Majorana levels, and (ii) much 
smaller than the gap Eg of the system, so that Majorana levels do not overlap in energy 
with levels of the bulk. In general, these two conditions depend on the system parameters 
in a nontrivial way given by the explicit form of 5m and Eg [6], namely, 

• 5m ^ e~^/^, where ^ is the localisation length of the Majorana-hke modes given by 

= min{|ln |a;+||, |ln|a;_||} with x± = ^^^2^(j+^\A\-r ^' 

• Eg ^2J - if |A| > J or if |A| < J and 2 J - |//| < 2|A|7J, and Eg = |A|(4 - 
l/ip/ (J^ — I Ap))^/^ otherwise (in the topologically nontrivial regime with 2J— > 0). 

Away from the hmit 2 J — — > where the gap closes (and where the topological phase 
transition occurs, in the limit N — >■ oo), (5m ^ T is generally more restrictive than Eg ^ F. 
As long as the above inequality J, |A| » F is ensured, both of these conditions can always 
be satisfied by tuning fi well inside the topologically nontrivial region, as shown in the 
numerical results presented in the main text. J, |A| ^ F is therefore a sufficient condition 
for our detection scheme to work, up to the tuning of /i. This condition is also necessary, 
since J, |A| ~ F would make the broadening comparable to the gap. 

Although this is not a requirement for our proposal, we finally remark that one typically 
expects to be able to work in a regime where J ~ |A|, since both quantities are essentially 
determined by the overlap of the Wannier modes of neighboring cavities. We also emphasize 
that the fabrication of a linear chain of ~ 10 almost identical cavities is technologically 
feasible, as suggested by the recent reahzation of a Kagome lattice of more than 200 nearly 
identical cavities [5, 7]. Photon correlation measurements on microwave cavities have also 
been recently demonstrated [8]. 



SIMULATING REAL-SPACE BRAIDING IN ID USING TUNNEL-BRAID OPER- 
ATIONS 



The ability to exchange — or "braid" — Majorana modes in real space plays an impor- 
tant role in demonstrating the non-Abelian nature of the latter. In the framework of the 
optical proposal presented in the main text, however, one is fundamentally restricted to one- 
dimensional systems where exchanging Majorana modes is physically impossible (Majorana 
modes would overlap and split in energy during the exchange, thereby hybridizing into com- 
plex fermionic modes and losing their exotic nature). It thus appears crucial to find another 
way to perform non-Abelian operations in the degenerate ground-state subspace associated 
with Majorana modes while preserving their degeneracy. Such a possibility was originally 
introduced in Ref. [9], where it was shown that so-called "tunnel-braid" operations can sim- 
ulate the braiding of Majorana modes in real space. In what follows, we translate these ideas 
into the framework of our optical proposal, and give an explicit procedure to simulate the 
exchange of two Majorana-like modes of light tunnel-coupled to the same ancillary cavity. 




L I R 



FIG. 1. Basic one-dimensional setup required to simulate non-Abelian operations such as braiding. 
Two Kitaev chains — left (L) and right (R) — with pump phases (j)L and respectively, exhibit 
Majorana-like modes 7^ and 7^ at their ends, as depicted. Each chain tunnel-couples with am- 
plitude Jl {Jr) to an intermediate nonlinear cavity (I) via its end cavity where 7^ (jr) is mostly 
localized. 

To illustrate the working principles of tunnel-braid operations, we consider a simple 
scenario where two Kitaev chains — left (L) and right (R) — are end-to-end tunnel-coupled 
through a single intermediate cavity (I) [Fig. 1]. We assume that the parametric pumps 
driving both chains have phases (pL and (ppt, respectively, and the same frequency Up. The 
intermediate cavity is nonlinear — though not driven — and resonant with Up, and both chains 
are in a topologically non-trivial phase with Majorana-like zero modes 7^ and 7^ located on 
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the left and on the right of the intermediate cavity, respectively. Photon tunneling occurs 
between the intermediate cavity and the end cavities of the Kitaev chains, as captured by 
the Hamiltonian 



Ht = -2JL{b[bi + h.c.) - 2jR{h\hR + /i.e.) , (21) 

where 6^,, hi and hR are annihilation operators associated with photons in the respective 
cavities (with obvious notations) , while Jl and Jr are positive amplitudes for tunneling into 
the left and right chains, respectively (factors of 2 are introduced for later convenience). 
Transforming to the fermionic picture defined in the main text, we obtain 

Ht = -2JL{a{ai + h.c.) - 2jR{a\aR + h.c), (22) 

with fermionic operators and aR that can be decomposed as 

a. = ie-^-/^(7i + i7L), ^^^^ 
aR^\e-"f-/\^R + \ij,), 

where 7^ and 7^ are Majorana operators whose explicit form is not relevant here. Since we 
are only interested in the low-energy physics associated with the Majorana-like zero modes 
7l and 7^, we can make in Eq. (22) the following replacement: 

(24) 

aR ^ |e-'*«/^i^, 
thus obtaining the low-energy effective tunnel Hamiltonian 

//t,eff = Jl ((-ie-'^-/^)4 - {-\e-"^-'yaj) 7l 

; / (25) 

+ Jr [{-e-"^^/^)a\ - {-e-"^^/yai) ^r. 

This can be recast in a more convenient form by performing the gauge transformation 
ai (-ie~''^^/^)a7, leading to 

/^t,eff = JM - a,)7L + JR{e--<f--^-+^y^a\ - e'^f--'^-+-y'aj)jR. (26) 

In this form, i^t.eff clearly appears asymmetric with respect to the interchange 7^ •<-> 7^? if 
the phase difference A0 = 0j? — 0l does not satisfy A0 = {An — l)7r for some integer n. 
Away from these fine-tuned values of A0, the tunnel Hamiltonian thus provides a means of 
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distinguishing between the different states of the Majorana qubit composed of 7^ and 7^, i.e. 
between different parity sectors of the full system. In order for this information to remain 
unaccessible to the system — thereby restoring what is usually referred to as the "topological 
protection" or "parity protection" of the Majorana qubit — it is therefore necessary to fine- 
tune the phase difference between the Majorana-like modes to some appropriate value (which 
can easily be done in our optical setting by adjusting the relative phase of the parametric 
pumps of each Kitaev chain; see main text). Assuming that this is satisfied — so that the 
even- and odd-parity sectors of the full system remain degenerate — ^the low-energy effective 
tunnel Hamiltonian reduces to 

= JLR{a\ - aihLR, (27) 

where Jlr — ^/ Jl + Jr and 'Ylr — {JlIl + JrIr) / Jlr- In this parity-protected scenario, 
everything therefore happens as if a single Majorana mode ^lr was tunnel-coupled to the 
intermediate cavity. Taking into account the possible detuning 5i of the intermediate cavity 
with respect to the frequency Up of the parametric pumps driving both Kitaev chains (see 
main text), the low-energy effective Hamiltonian describing the Majorana system tunnel- 
coupled to the intermediate cavity becomes 

i^eff 8ia\ai Jlr(4 - aihLR, (28) 

where, of course, \dj\, Jlr ^ Eg must be assumed (Eg being the energy gap of the Kitaev 
chains) in order for the low-energy description to be valid. It was shown in Ref. [9] that the 
above Hamiltonian allows for a variety of unitary "tunnel-braid" operations within the sub- 
space associated with the Majorana modes jl and 7^. In particular, starting from an empty 
intermediate cavity with Sj > and 5i/Jlr ^ 1, one can generate an adiabatic transition 
from an empty to an occupied intermediate cavity and a unitary rotation U{Jl, Jr) = j^r in 
the Majorana subspacc by adiabatically tuning 61/ Jlr from large positive to large negative 
values (keeping |5/|,./l_r, ^ Eg). The same unitary operation can be obtained if, starting 
from an occupied intermediate cavity with Sj < and \Si\/Jlr ^ 1, one adiabatically 
raises Sj/Jlr from large negative to large positive values. In both cases, the intermediate 
cavity and the Majorana system become entangled, and an inversion of the occupation of 
the intermediate cavity as well as a rotation ^lr in the Majorana subspace are obtained as 
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the latter return to an unentangled state at the end of the process. Using similar adiabatic 
processes, one can more generally perform non-Abelian operations within the degenerate 
ground-state subspace associated with 7^ and 7^. In particular, one can generate the op- 
eration Ub = + IlIr) and thus effectively simulate the braiding of 7^, and 7^ in real 
space (see, e.g., Ref. [10]) following the procedure below: 

1. Start with = Jr = (no tunnel coupling) and the intermediate cavity far blue- 
detuned {6 J ^ 0), such that it is initially empty. 

2. Switch on the tunnel coupling Jl, raising it to a constant value J <^ Sj. At this point, 
the intermediate cavity and the Majorana system are still effectively decoupled. 

3. Adiabatically red-shift the intermediate cavity, taking 5/Jlr from large positive to 
large negative values. In doing so, the intermediate cavity adiabatically changes from 
empty to occupied, and the unitary operation U{Jl = J, Jr = 0) = 7l is performed 
in the Majorana subspace. 

4. Switch on the tunnel coupling Jr, raising to the constant value Jl — J- 

5. Adiabatically blue-shift the intermediate cavity, taking 5/Jlr from large negative to 
large positive values. In doing so, the intermediate cavity adiabatically changes from 
occupied to empty, and the unitary operation U{Jl — J,Jr — J) — -^{il + 1r) is 
performed in the Majorana subspace. 

6. Switch off both Jl and Jr. 

At the end of the above cycle, the intermediate cavity and the Majorana system are left 
unentangled as they were initially. While the intermediate cavity returns to its initial empty 
state, a non-trivial unitary operation Ub — '^{11 + 1b)1l — IlIr) is obtained in 

the Majorana subspace, as if 7^ and 7^ had physically been exchanged. This exemplifies 
the fact that one can simulate real-space braiding in purely one-dimensional systems. We 
remark, however, that the tunnel-braid operations used in the above procedure are only 
parity-protected provided that the phase difference between the Majorana modes ^l and 
7ij satisfies A0 — (4n — l)7r (for some integer n). Because of this fine-tuning requirement, 
tunnel-braid operations are not strictly speaking topologically protected. Moreover, even if 
is fine-tuned so that the degeneracy of the Majorana modes is preserved, one must still 
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be able to control the tunneling amplitudes Jl and Jr sufficiently well in order to ensure that 
the desired tunnel-braid operations are generated. Indeed, different relative strengths Jl/ Jr 
generally lead to different operations (see Ref. [9] for further details). Although the lack 
of topological protection renders tunnel-braid operations somewhat unattractive for solid- 
state proposals where fermion parity conservation is almost perfect, it does not constitute 
an additional problem in our optical setting where parity is anyway only conserved on 
timescales much shorter than the lifetime of the photons. If tunnel-braid operations can be 
completed within such timescales, the non-Abelian nature of Majorana-like modes of light 
can be simulated. 
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